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Abstract 

The best motivated alternatives to general relativity are scalar-tensor theories, in which the 
gravitational interaction is mediated by one or several scalar fields together with the usual gravi- 
ton. The analysis of their various experimental constraints allows us to understand better which 
features of the models have actually been tested, and to suggest new observations able to discrim- 
inate between them. This talk reviews three classes of constraints on such theories, which are 
qualitatively different from each other: (i) solar-system experiments; (ii) binary-pulsar tests and 
future detections of gravitational waves from inspiralling binaries; (hi) cosmological observations. 
While classes (i) and (ii) impose precise bounds respectively on the first and second derivatives 
of the matter-scalar coupling function, (iii) a priori allows us to reconstruct the full shapes of 
the functions of the scalar field defining the theory, but obviously with more uncertainties and/or 
more theoretical hypotheses needed. Simple arguments such as the absence of ghosts (to guarantee 
the stability of the field theory) nevertheless suffice to rule out a wide class of scalar-tensor mod- 
els. Some of them can be probed only if one takes simultaneously into account solar-system and 
cosmological observations. 



* Invited talk at the workshop Phi in the Sky: The Quest for Cosmological Scalar Fields, Porto, July 2004 
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I. INTRODUCTION 



In the most natural alternative theories to general relativity (GR), gravity is mediated not 
only by a (spin-2) graviton corresponding to a metric g^, but also by a (spin-0) scalar field 
(p. Such scalar partners generically arise in all extra-dimensional theories, and notably in 
string theory. A dilaton is indeed already present in the super mult iplet of the 10-dimensional 
graviton, and several other scalar fields (called the moduli) also appear when performing 
a Kaluza-Klein dimensional reduction to our usual spacetime. They correspond to the 
components of the metric tensor g mn in which m and n label extra dimensions. Moreover, 
contrary to other alternative theories of gravity, scalar-tensor theories respect most of GR's 
symmetries: conservation laws, constancy of non-gravitational constants, and local Lorentz 
invariance even if a subsystem is influenced by external masses. They can also satisfy exactly 
the weak equivalence principle (universality of free fall of laboratory-size objects) even for a 
massless scalar field. 

Scalar fields are also involved in the cosmological models which reproduce most consis- 
tently present observational data. In particular, inflation theory is based on the presence of 
a scalar *p in a potential V((p) (for instance parabolic). It behaves as a fluid with a positive 
energy density SnGp^, = (p 2 + 2V(<p) but a negative pressure SnGp^ = (p 2 — 2V(<p). This 
causes a period of exponential expansion of the universe, which can explain why causally 
disconnected regions at present may have been connected long ago. The isotropy of the ob- 
served cosmic microwave background (CMB) can thus be understood. Inflation also predicts 
that our universe is almost spatially flat, because any initial curvature has been exponen- 
tially reduced by the expansion. This is in remarkable agreement with the location of the 
first acoustic peak of the CMB spectrum at a multipolar index I ~ 220 [l| . Various observa- 
tions, notably of type la supernovae |2j, tell us that there is about 70% of negative-pressure 
dark energy in our present universe (Q\ ~ 0.7), suggesting that its expansion has been 
re-accelerating recently, since redshifts z ~ 1. This can be explained by the presence of a 
cosmological constant A in GR, but the quantity Vt\ ~ 0.7 translated in natural units gives 
an extremely small value A ~ 3 x 10~ 122 c 3 /(hG), very problematic for particle physics if 
A is to be interpreted as the vacuum energy. This is the main reason why "quintessence" 
models have been proposed, in which the cosmological constant is replaced again by the po- 
tential V(ip) of a scalar field. Its evolution towards a minimum of V during the cosmological 
expansion then explains more naturally why the present value V^(^o) — A/2 is so small. 

Besides these theoretical and experimental reasons for studying scalar-tensor theories of 
gravity, one of their greatest interests is to embed GR within a class of mathematically 
consistent alternatives, in order to understand better which theoretical features have been 
experimentally tested, and which can be tested further. To simplify this review, we will 
restrict it to models involving a single scalar field, although the study of tensor-multi-scalar 
theories can also be done in great detail [3j]. It suffices to note that their phenomenology 
is richer but similar to the single scalar case, at least when all scalar degrees of freedom 
carry positive energy, as required for the vacuum to be stable. In view of the very precise 
experimental verifications of the weak equivalence principle, we will also focus this review on 
theories which satisfy it exactly. In other words, we will assume that the action of geometry 
on matter is the same as in GR, but that the dynamics of geometry and the action of matter 
on it is modified because of the presence of a scalar field (see J. P. Uzan's contribution to the 
present proceedings, and in particular the upper-right panel of its Fig. 1). All matter fields, 
including gauge bosons, will thus be assumed to be universally coupled to one second rank 
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symmetric tensor, say g^ u . The difference with GR lies in the kinetic term of this tensor, 
which is now a mixing of spin-2 and spin-0 excitations. One indeed defines g^ v = A 2 ((p)g fJiU , 
where A(ip) is a function of the scalar field, while g^ v denotes the Einstein (spin-2) metric. 
We will mainly consider the following action 



which depends on two functions of the scalar field: the coupling function A(ip) to matter, 
and a potential V(ip). [Our signature is -+++, R is the scalar curvature of g^, and g its 
determinant.] 

In Section [HI we will review the main experimental constraints on this coupling function 
A(<p), coming from two qualitatively different classes of data: solar-system experiments 
and binary-pulsar observations 0, 0- in Section ITTTl we will underline that cosmological 
observations give again very different types of constraints: Although they are more noisy, 
they a priori allow the reconstruction of both A(cp) and V((p) U E|. The most general 
scalar-tensor theories will be defined in Sec. II VI and we will illustrate in Sec. 13 on a very 
special but interesting case, that one must sometimes take simultaneously into account solar- 
system and cosmological observations to obtain significant constraints on a theory |9[. We 
will finally give our conclusions in Sec. IVII 

II. SOLAR-SYSTEM AND BINARY-PULSAR CONSTRAINTS 
A. Solar-system tests 

A massive scalar field has a negligible effect on the motion of celestial bodies if its mass is 
large with respect to the inverse of the interbody distances. On the other hand, if its mass 
is small enough, its potential V(ip) can be locally neglected, but its coupling function to 
matter, A(<p), is strongly constrained by experiment, as we will see below. The intermediate 
case, for which the size of the experiment is of the same order as the inverse scalar mass, 
is the main topic of E. Adelberger's contribution to the present proceedings (see also (Io|). 
The "chameleon" field [llj is another case where both the potential and the matter-scalar 
coupling function must be taken into account in solar-system experiments; see P. Brax's 
contribution to these proceedings. In the present section, we will review the constraints on 
A(ip) when the potential V(ip) can be neglected. 

The precision of solar-system observations allow us not only to test Newton's law, but also 
its relativistic corrections oc 1/c 2 , called "post-Newtonian". At this order, the predictions of 
metric theories of gravity can be parametrized by a set of 10 real numbers in the so called 
"PPN" formalism (parametrized post-Newtonian) [12j. All of them are presently constrained 
to be very close to their general relativistic values (l3|, and in particular the two famous 
parameters /3 PPN and 7 PPN introduced by Eddington in the Schwarzschild metric: 
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FIG. 1: Solar-system and binary-pulsar constraints on the matter-scalar coupling function 
In A(ip) = ao((f — (fo) + ^Po((f — fo) 2 + 0(f — <po) 3 - The allowed region is shaded. The ver- 
tical axis (/?o = 0) corresponds to Brans-Dicke theory with a parameter 2wbd + 3 = 1/ag- The 
horizontal axis («o = 0) corresponds to theories which are perturbatively equivalent to GR, i.e., 
which predict strictly no deviation from it (at any order l/c n ) in the weak- field conditions of the 
solar system. 



General relativity corresponds to /3 PPN = 7 PPN = 1. In scalar-tensor theories 0,0, these 
two latter parameters are the only ones which can differ from their GR values. They are 
related to the first two derivatives of 

In A((p) = a & - (p ) + ^[3 ((p - <^ ) 2 + C% - <^ ) 3 , (3) 

computed at the background value (p of the scalar field: 



I'l'x , 2a o ^ppn -, 1 aoPoap 



The factor ckq comes from the exchange of a scalar particle between two bodies, whereas 
«oA)tto comes from a scalar exchange between three bodies. Solar system experiments 
impose the constraints displayed as thin lines in Fig. ^ The Mercury symbol refers to the 
perihelion shift of this planet, whose observed value imply the bound [14J 

|2 7 PPN_ /? PPN_ 1 | <3xl0 -3 ; (5) 

the Moon symbol refers to Lunar Laser Ranging p~5| : 

4/? PPN _ 7 PPN _ 3 = ( _ Q ? ±1)y< 10 - 3) (6) 
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FIG. 2: Heuristic argument to explain the phenomenon of "spontaneous scalarization" . When 
Po < and the compactness Gm/Rc 2 of a body is large enough, it is energetically favorable to 
create a local scalar field different from the background value. The body becomes thus strongly 
coupled to the scalar field. 

the star symbol to light deflection as measured by Very Long Baseline Interferometry [lrjl ]: 

| 7 PPN - 1| <4x 1(T 4 , (7) 

and the label "Cassini" to the impressive recent constraint obtained by measuring the time 
delay variation to the Cassini spacecraft near solar conjunction |l7| |: 

7 ppn_x = (2.1 ±2.3) x irr 5 . (8) 

Solar-system tests thus constrain the first derivative «o to be small, i.e., the linear interaction 
between matter and the scalar field to be weak. On the other hand, the second derivative 
Po may take large values (of either sign), i.e., matter may be strongly coupled to two scalar 
lines. 



B. Nonperturbative strong- field effects 

At higher post-Newtonian orders l/c n , a simple diagrammatic argument shows that any 
deviation from GR involves at least two factors ao, and has the schematic form 

deviation from GR = «qX 

where m and R denote the mass and radius of the considered body, and Ao, Ai, ... are 
constants built from the coefficients ao, Po, ... of expansion (J3J). Since is experimentally 
known to be small, we thus expect the theory to be close to GR at any order. However, some 
nonperturbative effects may occur in strong-field conditions: If the compactness Gm/Rc 2 
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FIG. 3: Scalar charge «a versus baryonic mass m\, for the model A(tp) = exp(— 3<p 2 ) (i.e., (3q = 
—6). The solid line corresponds to the maximum value of «o allowed by Eqs. and the 

dashed line to «o = 0- The dotted lines correspond to unstable configurations of the star. 



of a body is greater than a critical value, the square brackets of Eq. can become large 
enough to compensate even a vanishingly small Oq. This can happen notably for neutron 
stars, whose compactnesses are of order Gm/Rc 2 ~ 0.2, as compared to only 2 x 10 -6 for 
the Sun or even 7 x 10 -10 for the Earth. To illustrate this, let us consider a model for 
which ao vanishes strictly, i.e., which is perturbatively equivalent to GR: There is strictly no 
deviation from GR at any order in a perturbative expansion in powers of 1/c. A parabolic 
coupling function ha.A(ip) = |/3 </9 2 suffices for our purpose (we set here <p = to simplify). 
At the center of a static body, the scalar field takes a particular value ip c , and it decreases 
as 1/r outside. The energy of such a scalar field configuration involves two contributions, 
coming respectively from the kinetic term and from the matter-scalar coupling function in 
action (JXJ) . As a rough estimate of its value, one can write 



Energy 



™ 2 (4^ +eft * 2) - (io) 



When (3q < 0, this is the sum of a parabola and a Gaussian, and if the compactness Gm/Rc 2 
is large enough, the function Energy (ip c ) has the shape of a Mexican hat, see Fig. |2j The 
value ip c = now corresponds to a local maximum of the energy. It is therefore energetically 
favorable for the star to create a nonvanishing scalar field (f c , and thereby a nonvanishing 
"scalar charge" d\n A({p c )/d<p c = j3o<p c - This phenomenon is analogous to the spontaneous 
magnetization of ferromagnets. 

This heuristic argument has been verified by explicit numerical calculations, taking into 
account the coupled differential equations of the metric and the scalar field, and using various 
realistic equations of state to describe nuclear matter inside a neutron star Q . The correct 
definition of the linear coupling strength between a compact body A and the scalar field 
reads a a = dlnm^/dipo. It is plotted in Fig. El for the particular model (3$ = —6. One finds 
that there exists indeed a "spontaneous scalarization" above a critical mass (whose value 
decreases as — (3 grows). 

On the other hand, if (3q > 0, both the above heuristic argument and the actual numerical 
calculations show that |oa| < |«o|- I n that case, one finds that neutron stars are even less 
coupled to the scalar field than solar-system bodies. This decoupling is similar to the 



behavior of the "chameleon" field but it is here caused by the nonlinearity of lnA(<^), 
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FIG. 4: Mass plane (tua = pulsar, me = companion) of the Hulse- Taylor binary pulsar PSR 
B1913+16 in general relativity (left panel) and for a scalar-tensor theory with /3q = —6 (right 
panel). The widths of the lines are larger than la error bars. While GR passes the test with flying 
colors, the value (5q = —6 is ruled out. 

whereas the chameleon's decoupling is due to a competition between A(ip) and the potential 
V(tp). 

C. Binary-pulsar tests 

The scalar charge a a enters the predictions of the theory in the same way as «o m weak- 
field conditions. For instance, in the orbital motion of two bodies A and B, the Eddington 
parameter -)/ PPN keeps the form of Eq. but it now involves the product «a«b of the two 
scalar charges instead of aft. Similarly, the strong- field analogue of /3 PPN involves products of 
scalar charges and their derivatives /3a = da\/d<p . Since |«a| — 0.6 in the model of Fig. 01 
one thus expects deviations by ~ 35% from some general relativistic predictions. Moreover, 
the quadratic coupling strength (3a can take very large numerical values near the critical 
mass, like the magnetic susceptibility of ferromagnets. Therefore, even larger deviations 
from GR are found when the mass of a neutron star happens to be close to the critical one. 

Binary pulsars are ideal tools for testing gravity in strong-field conditions. A pulsar is a 
rapidly rotating neutron star emitting a beam of radio waves, like a lighthouse. Experiment 
tells us that isolated pulsars are very stable clocks, when they are old enough. A pulsar 
A orbiting a companion B is thus a moving clock, the best tool that one could dream 
of to test a relativistic theory. Indeed, by precisely timing its pulse arrivals, one gets a 
stroboscopic information on its orbit, and one can measure several relativistic effects. Such 
effects do depend on the two masses tua, tobj which are not directly measurable. However, 
two different effects suffice to determine them, and a third relativistic observable then gives 
a test of the theory. In the case of the famous Hulse- Taylor binary pulsar PSR B1913+16, 
three relativistic parameters have been determined with great accuracy 0] : (i) the Einstein 
time delay parameter jt, which combines the second-order Doppler effect (oc v\/2c 2 , where 
va is the pulsar's velocity) together with the redshift due to the companion (oc Gitis/tabc 2 , 
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where tab the pulsar-companion distance); (ii) the periastron advance u> (relativistic effect 
of order v 2 /c 2 ); and (iii) the rate of change of the orbital period, P, caused by gravitational 
radiation damping (an effect of order v 5 / c 5 in GR, but of order v 3 / c 3 in scalar-tensor theories; 
see below). Figure |U displays the plane of the two a priori unknown masses tua and m-Q. For 
each relativistic parameter, the prediction of a given theory is consistent with the observed 
value only along a thin line. In GR, the fact that the three lines meet in one point means that 
there exists a pair of masses {m^m-o) simultaneously consistent with all three observables. 
This is thus a spectacular confirmation of Einstein's theory. 

Obviously, these lines are deformed in scalar-tensor theories, and if they do not have 
any common intersection, the corresponding model is ruled out. The right panel of Fig. 0] 
illustrates such a case. The allowed theories lie below and to the right of the line labeled PSR 
B1913+16 in Fig. ^ This plot shows vividly the qualitative difference between solar-system 
and binary-pulsar observations. Indeed, the latter impose 

/?o>-4.5, (11) 

even for a vanishingly small cuo- This constraint is due to the spontaneous scalarization of 
neutron stars, which occurs when — (3q is large enough. Equations allow us to rewrite this 
inequality in terms of the Eddington parameters /3 PPN and 7 PPN , which are both consistent 
with 1 in the solar system. One finds 

aPPN _ 1 

^PN3I< 1 - 1 - ( 12 ) 

The singular (0/0) nature of this ratio underlines why such a conclusion could not be ob- 
tained in weak-field experiments. 

Several other binary pulsars are presently known and accurately timed. To simplify Fig.Q 
we did not display all the constraints they impose in the theory plane (see Ref. p). We just 
plotted as a dotted line the ti ght ones imposed by the recently timed neutron star-white 
dwarf binary PSR J1141— 6545 [19]. Contrary to neutron stars, the white dwarf companion is 
never spontaneously scalarized, because its compactness Gm/Rc 2 is too small. This binary 
system is thus very asymmetrical, and the pulsar's scalar charge is generically much 
greater than its companion's, «b- This causes a huge deviation from GR in the orbital 
period variation P. Indeed, the energy flux carried out by gravitational waves is of the form 
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The first curly brackets contain the prediction of general relativity, of order v 5 /c 5 , whereas 
the second ones contain the extra contributions predicted in tensor-scalar theories. In par- 
ticular, the dipolar contribution is of order v 3 /c 3 , much larger than the usual quadrupole 
of GR if the two scalar charge are significantly different. This is the reason why the asym- 
metrical system PSR J1141— 6545 is so constraining for scalar-tensor theories, although its 
experimental uncertainties on P are still rather large. Notice that this binary pulsar is al- 
most as constraining as solar-system tests even in the region /3q > 0; see Fig. ^ It should 



probe values of the Eddington parameter |7 PPN — 1| ~ 10~ 6 by the end of the decade, i.e., 
an order of magnitude better than present solar-system limits (|5|). 

The LIGO and VIRGO interferometers will detect gravitational waves emitted by inspi- 
ralling binaries. Since matter is known to be weakly coupled to the scalar field in the solar 
system (small value of |a |), these detectors will not be sensitive of the helicity-0 waves. 
On the other hand, the time evolution of the helicity-2 chirp does depend on the energy 
flux (|13|) . which differs significantly from the GR prediction when the scalar charges a a, 
ctB are nonzero. Therefore, the GR wave templates used for matched filtering in LIGO 
and VIRGO may not be accurate if there exists a scalar partner to the graviton, and the 
signal-to-noise ratio may then drop. Fortunately, it was shown in |5( that binary-pulsar data 
are so precise that they already exclude the models which would have predicted significant 
effects in the gravitational waveforms. Therefore, although these interferometers are a priori 
more sensitive to the scalar field than classic solar-system tests, one may securely use the 
GR wave templates for their data analysis. On the other hand, it was shown in |2(| that 
the LISA space interferometer can be sensitive to scalar effects which are still allowed by all 
present tests. Future binary-pulsar data should nevertheless probe them before the LISA 
mission is launched p. 

In conclusion, solar-system tests tightly constrain the first derivative of In A(ip) (linear 
matter-scalar coupling strength a ), whereas binary-pulsar data impose that its second 
derivative (5q (quadratic coupling matter-scalar-scalar) is not large and negative. We will now 
see that cosmological observations give access to the full shape of this coupling function, of 
course not with the same accuracy as the above tests, but with the capability of constraining 
any higher derivative of In A(tp) (vertex of matter with any number of scalar lines) . Moreover, 
cosmological data can also give access to the full shape of the potential V((p). 



III. RECONSTRUCTION OF A SCALAR-TENSOR THEORY FROM COSMO- 
LOGICAL OBSERVATIONS 

In cosmology, the usual approach to study quintessence models is to assume a particular 
form for the potential V(ip) (and the matter-scalar coupling function A(ip) when one con- 
siders "extended quintessence" models), to compute all possible observable predictions, and 
to compare them to experimental data. 

In contrast, in the phenomenological approach, one wishes to reconstruct the Lagrangian 
of the theory from cosmological observations. It was proved in [7( that the knowledge of the 
luminosity distance Dl(z) and of the density fluctuations S m (z) = 5p/p as functions of the 
redshift z indeed suffices to reconstruct both the potential V(ip) and the coupling function 
A((p). [The knowledge of the Hubble constant H or the present matter density fi m o is not 
even necessary, as they derive from Dl(z) and 5 m (z).] In order to write the corresponding 
equations as simply as possible, it is convenient to perform a change of variables on action (JJ). 
Let us define 

g ilv = A 2 (ip)g tlv ; ^ = A-\ip); 
2uj BD ^) + 3 = A 2 (cp)/A' 2 ^); U($) = 2V(<p)/A 4 (<p) . (14) 

Then action (Q) takes the "Brans-Dicke" form 

S = T^ \ ^V^(^-^^(^$) 2 -2^ (15) 
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in which the scalar curvature R, the determinant 7j, and the contraction of indices in (c^$) 2 
now all correspond to the physical ( "Jordan- frame" ) metric . Note that this action defines 
exactly the same theory as Eq. (ffl). 

The reconstruction equations can then be written straightforwardly. The first step is to 
derive the Hubble function H(z) from the knowledge of the luminosity distance Dl(z)\ 
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where a prime denotes here a derivation with respect to the redshift z, and Q Kt o is the possible 
(small) contribution of the Universe's spatial curvature to the total energy density. In the 
short wavelength limit A <C {H~ l , m^ 1 ), the field equations for small perturbations yield 

5 m + 2H5 m — 4:7iG e spS m ~ 0, where G c s = G[A 2 (ip) + (dA/dip) 2 } is the effective gravitational 
constant between two close test masses, taking into account the scalar-field contribution. In 
terms of the Brans-Dicke variables, this gives 
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On the other hand, the background evolution equations yield 
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Therefore, Eqs. (|T7 |l -(fIU jl allow us to reconstruct ®(z), U(z) and u B d(z) as functions of the 
redshift, and thereby U($>) and ubd{&) in a parametric way. Actually, the last factor of 
Eq. (fT7j) may be neglected, since we know that (2u;bd + 3) _1 = a 2 , 1 at present. In that 
case, Eqs. (|17j) - (jl9|) become merely algebraic, and no integration of differential equation is 
required. However, one should note that second derivatives of experimental data are needed, 
so that large uncertainties are expected in such a reconstruction program (see R. Jimenez's 
contribution to the present proceedings). The point that we wish to emphasize is that the 
full reconstruction of the microscopic action ()15|) is theoretically possible. Although it needs 
some algebra, this result seems anyway obvious: It is possible to fit two observed functions 
[D L (z) and S m (z)] thanks to two unknown ones \V(ip) and A(ip), or equivalently U(§) and 

«>bd($)]- 

However, future experiments (like the SNAP satellite) will only give access to the lu- 
minosity distance D L (z) with a good accuracy, and the density contrast S m (z) cannot yet 
be used to constrain the models. A semi-phenomenological approach can thus be useful: 
We make some theoretical hypotheses on either the potential V((p) or the coupling function 
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FIG. 5: Matter-scalar coupling function In A((p) which reproduces exactly the same accelerated 
expansion as GR plus a cosmological constant (with Q\ ~ 0.7), without any scalar field potential 
V(<p) nor any cosmological constant. However, this model works only for redshifts z < 0.7. 

A(<p), and we reconstruct the other one from Dl(z). A priori, one may think that such a 
reconstruction is again obvious: We fit one observed function with one unknown 

function (V(<£>) or But this naive reasoning is only valid locally, on a small interval. 

Indeed, the reconstructed function may for instance diverge for some value of the redshift, 
or one of the degrees of freedom may need to take a negative energy beyond a given redshift, 
which would make the theory unstable (and ill defined as a field theory on the surface where 
the energy changes its sign). The positivity of the graviton energy implies A 2 ((p) > 0, which 
can be translated in terms of the Brans-Dicke scalar field as $ > 0. On the other hand, the 
positivity of the scalar-field energy imposes the minus sign in front of the scalar kinetic term 
— {d^ip) 2 in action which translates as c^bd > — § m terms of the Brans-Dicke parameter. 
It was shown in Ref. 0] that these conditions impose tight constraints on the theories as 
soon as one knows Di{z) over a wide enough interval z G [0, ~ 2]. 

For instance, the present accelerated expansion of the universe can be perfectly described 
by a scalar-tensor theory with a vanishing potential V((p) = (and therefore a vanishing 
cosmological constant too). Reference [H derived analytically the coupling function A(<p) 
which reproduces exactly the same evolution of the scale factor a(z) as the one predicted 
by GR plus a cosmological constant. Figure displays this reconstructed function \nA(<p), 
which has a nice parabolic shape, with a minimum very close to the present value <pq of the 
scalar field, and a positive second derivative. This is not only consistent with binary-pulsar 
data [which forbid large and negative values of this second derivative /3 , see Eq. (Jll|) ] but 
also with the cosmological attractor phenomenon analyzed by Damour and Nordtvedt |2l| : 
The scalar field is generically attracted towards a minimum of In A(ip) during the cosmolog- 
ical expansion, whereas some fine tuning would be necessary to reach a maximum (negative 
/3q). Therefore, we are in the difficult situation in which two very different theories are both 
consistent with experimental data, and there seems to be no way to distinguish them. For- 
tunately, Ref. found that this scalar-tensor theory cannot mimic GR plus a cosmological 
constant beyond a redshift z ~ 0.7, because the scalar field cp would diverge at this value, 
and above all because the graviton energy would become negative beyond. Therefore, it 
suffices to measure D L (z) precisely enough up to z ~ 1 to rule out such a potential-free 
scalar-tensor theory. Actually, if Dl(z) is measured over a wider interval z G [0, ~ 2], large 
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experimental errors (tens of percent) are not problematic: It is still possible to distinguish 
this potential- free model from GR plus a cosmological constant, and thereby to rule out 
one of them. The results of the SNAP satellite up to z ~ 2 will therefore be very useful to 
constrain scalar-tensor theories of gravity. 

One can also impose a particular form of the coupling function A(ip) and reconstruct the 
potential V((p) which reproduces the observed luminosity distance Dl(z). For instance, for 
a minimally coupled scalar field A(ip) = 1 (usual "quintessence") in a spatially curved uni- 
verse, Ref. analytically derived the expression of V(cp) which gives the same cosmological 
evolution as GR plus a cosmological constant in a spatially flat universe. It was found that 
the shape of the potential is smoother when the universe is (marginally) closed. If it is flat or 
almost flat, one obviously recovers a cosmological constant with its unnaturally small value 
A~3x 10~ 122 c 3 /(hG). Therefore, in that case, aesthetic reasons may help us discriminate 
between the theories, instead of the much stronger argument of the positivity of energy that 
was used above. This shows anyway that the sole knowledge of Dl{z) suffices to constrain 
scalar-tensor theories of gravity. 



IV. THE MOST GENERAL SCALAR-TENSOR THEORIES 

More general scalar-tensor theories than action (JTJ) can be defined. First, one may con- 
sider several scalar fields 0], labeled by latin indices, (p a . In the action, their kinetic term 
— ^df J ,(p a d fl <p b will thus involve a matrix j a b to contract the indices. Such a matrix may itself 
depend on the scalar field, and it defines thus what is called a a-model. An example of a 
tensor-bi-scalar theory is provided by hybrid inflation. Note that the cr-model metric 
must be positive definite for all the scalar fields to carry positive energy. 

The quintessence potential V(ip a ) may also depend on the various scalar fields, and 
thereby define n-dimensional mountains along which the scalar degrees of freedom roll down. 
Similarly, the matter-scalar coupling function A((p a ) may have a nontrivial shape, and define 
a rich "extended quintessence" phenomenology. 

If the weak equivalence principle is not exactly satisfied, different species of matter, 
labeled i, may have different coupling functions Aj(y? a ) to the scalar field. It should be noted 
that superstring theory does predict such couplings, which violate equivalence principle 
tests by 9 or 10 orders of magnitude at tree level (see J. P. Uzan's contribution to the 
present proceedings). To reconcile it with experiment, the standard argument is that no 
symmetry protects the masses of the scalar fields, so that all of them should acquire a large 
enough mass at present to be exponentially suppressed, or anyway a potential such that 
they are very weakly coupled to matter on Earth (cf. the chameleon model [11]). However, 
"flat directions" generically arise in low-energy effective models inspired by superstrings, 
i.e., one or several scalar fields which remain strictly massless and violate a priori the 
weak equivalence principle. In that case, string loop corrections can induce an attractor 
mechanism 25l . which efficiently drives such scalars towards an extremum of their coupling 



functions to matter, so that they almost decouple at present. 

It is also possible to couple the scalar fields to the Gauss-Bonnet topological invariant, 
R 2 —4R 2 +R 2 . This is the only combination involving powers of the Riemann and/or the 
Ricci tensors which preserves the spectrum of the theory, i.e., which does not add any extra 
degree of freedom. A function f(R) of the scalar curvature is well known to be equivalent 
to an extra scalar degree of freedom in the model, cf. 2^1 . It suffices to rewrite f f(R) as 



J [f(ip) + (R — <f)f'(f)] and to show that the field equations deriving from these actions 
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coincide. The latter is of the Brans-Dicke type (JHJ), with ubd = 0, and a mere change of 
variables (conformal transformation) g^ u = f'{<p)g in , yields the canonical form ((TJ). Similarly, 
a function f(R, OR, . . . , n n R) of the scalar curvature and its iterated d'Alembertian can be 
shown to be equivalent to the introduction of n + 1 scalar fields in the theory [23|. On 
the other hand, functions JiR^u) or f{R^ vpa ) of the Ricci or the Riemann tensor have 
been proved several times [24] to involve generically an extra massive spin-2 ghost, i.e., an 
extra negative-energy graviton which spoils the stability of the theory. The only allowed 
combination is the Gauss-Bonnet topological invariant. 

Finally, one may also consider a function F^d^tpd^tp) of the kinetic terms, as in "k- 
essence" models |26|. However, only some particular functions F preserve the positivity of 
the scalars' kinetic energy. 

In conclusion, the most general tensor-multi-scalar theories that we may consider are 
defined by the following action: 

-hf V=9 W(<f) [Rl vpu - ARl + R 2 ) 

+ [matter, ; = A\{y a )g^\ . (20) 



V. EXPERIMENTAL CONSTRAINTS ON A SCALAR GAUSS-BONNET COU- 
PLING 

In this section, we will consider the very particular case of a scalar-tensor theory with a 
single scalar field ip, no matter-scalar coupling [A(ip) = 1], no potential [V(<p) = 0], and a 
standard kinetic term — |((9 M ^) 2 , but with a scalar-Gauss-Bonnet coupling W((p): 

s = ^gJ ^\^-\ (5^) 2 } - n J V^w( v ) {R\ upa - m% + r 2 ) 

+5'matter [matter ; g pu ] . (21) 

Our aim is just to illustrate the kind of constraints which can be imposed on such a scalar- 
Gauss-Bonnet coupling This model will underline that solar-system and cosmological 
observations must be sometimes simultaneously taken into account to constrain a theory. 
The perturbative but highly nonlinear effects which happen can also be instructive. 

Solar-system and binary-pulsar tests are local, and any deviation from GR depends on 
the magnitude of the scalar field created by a massive body. Let us thus analyze first the 
equation satisfied by <p in the vicinity of a spherical mass m®. We can assume that the 
metric is close to the Schwarzschild solution, and we get at the first nonvanishing order in 
powers of Gm Q / c 2 

^ = i^p) 2 W + W ^ + °^ 2 )} , (22) 

where we have set £q = 16irGh/c 3 , and where the derivative W'(ip) has been expanded in 
powers of <p in the right-hand side. Since we are assuming that p> takes small values, let 
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FIG. 6: Scalar-Gauss-Bonnet coupling function W(ip) which exactly reproduces the cosmological 
expansion predicted by GR plus a cosmological constant. 



us neglect the contribution W^'f. We can then compute any observable prediction, but we 
quote below only the results for the light deflection angle (A0*) and for the perihelion shift 
per orbit (A8 p ), which suffice for our purpose: 



= ^ + ^ ^? - K, (23) 



4Gm 1536 / Gm Q \ i c \ TXr/2 
r Q c 2 35 V r oC 2 / V r o, 

^=^ + l*,(^)'g)V (24, 

where tq is the minimal distance between the light ray and the Sun, and p is the semilatus 
rectum of an orbit. The first terms on the right-hand sides are the usual GR predictions, at 
first order in Gm Q /c 2 . In conclusion, solar-system (and binary-pulsar) tests can easily be 
passed if \Wq\ is small enough. 

One can now reconstruct the full shape of W(<p) from the cosmological observation of the 
luminosity distance Dl(z), as in Sec. 11111 above. We find that this can always been done, 
without any problem of negative energy, contrary to what happened above for the matter- 
scalar coupling function A((p). Moreover, there exists again an attraction mechanism which 
drives the scalar field towards a minimum of W(if) during the cosmological expansion. 
Therefore, a small value of the slope \Wq\ is indeed expected at present, consistently with 
what is needed for solar-system tests. In conclusion, we are faced again with a serious 
problem: We just found a theory which seems to be consistent with all experimental data, 
although it is very different from GR in its field content. Our aim is therefore to find a way 
to distinguish it from GR, or to rule it out for internal consistency reasons. 

Figure El displays this reconstructed coupling function W(ip), in which the present value 
of the scalar field is close to the minimum at ip = 0. Its shape is nicely smooth, but its 
second derivative at the origin is huge if one divides it by the tiny factor (167rG h/c b )Hl. 
One gets Wq ~7x 10 119 , which is in fact not surprizing, since the coupling function W(<p) 
behaves in action (|2*Tj) as the inverse of a cosmological constant. Indeed, W(<p) multiplies 
the square of the curvature tensor, whereas the usual Einstein-Hilbert term involves the first 
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power of the scalar curvature, and a cosmological constant does not multiply any curvature 
term at all. It was thus expected that W(<p) involve a dimensionless number of the order of 
the inverse of (HG/c 3 )A ~3x 10~ 122 . Therefore, this model is ugly, but it is not yet ruled 
out. One should not confuse fine tuning and large (or small) dimensionless numbers in a 
model. We are here in the second situation, but there is a priori no fine tuning since the 
scalar field is attracted towards the minimum of W(<p) during the cosmological expansion. 
There remains to study how efficiently it is attracted, but this is actually not necessary for 
our purpose. 

Indeed, Wq takes such a gigantic value that an approximation that we made to analyze 
solar-system tests is no longer valid. Indeed, we have |WqVI ^ I^oIj so that the second 
term on the right-hand side of Eq. (J22j) cannot be neglected. To simplify the discussion, we 
will anyway assume that W((p) is parabolic, which is a good approximation in a vicinity of 
the minimum (p — 0. We will thus neglect the higher order terms 0((p 2 ) in Eq. (J22j) : taking 
them into account would not change our conclusions below. We did not find a close analytic 
solution to Eq. (J22j) . but is is possible to write it as a series 

w = ]^ V 1 ( 12%&m%Wg \ n 

Wg^ (3 x4)(7x 8)---(4n- l)(4n) V ^ ) 1 ' 

(26) 



Wo 



cos / Gm Q £ 



cosh \ r 2 c 2 V^l^o 



3|WX'| - 1 



if WS < 0, 
if WS > 0. 



The second expression is a good approximation if the argument of the cosine (or hyperbolic 
cosine if WS > 0) is much greater than 1. This is the case if we use the huge value of Wq 
obtained above from the cosmological reconstruction, and a typical solar-system distance 
for the radius r: The argument of the hyperbolic cosine is then of order 10 s . 

The above solution is such that ip oc Wq, therefore we do not find any nonperturbative 
effect similar to the "spontaneous scalarization" of neutron stars mentioned in Sec. Ill Bl 
above. Moreover, ip — > as r — > oo, and we recover GR for distances r > 4 x 10 14 m 
(i.e., farther that the solar system including Oort's comet cloud). On the other hand, there 
are highly nonlinear corrections proportional to l/r 4ri within the solar system. Since the 
ratio (12£ qG 2 m^W g /r^c 1 ) is much greater that 1, its successive powers blow up, but they 
are compensated by the factors 1/(3 x 4 x 7 x ■• -4n) which behave like the inverse of 
factorials. Therefore, the successive terms of series (|23j) start to grow exponentially, then 
reach a maximum for a value of the index n which may be large, and finally tend towards 
zero. Each of these successive terms must be assumed to be small enough for the model to 
pass all classical tests, but one should not forget that the largest one does not correspond 
to n — 1. 

In order to study the effects of such highly nonlinear terms in the solar system, we 
compute their corrections to the Schwarzschild metric in the form 



fi + E§) cV + f 1 + E dr2 + r2 i d ° 2 + sin2 6 > ( 27 ) 

\ n>l J \ n>l / 



and we find that the light deflection angle and the perihelion shift are respectively given by 



r (n + 1 

n>l v 
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QnGm Q ^ n(n-l)f3 n c 2 2 
Mp = ~ ^ 2Gm eP ^ * + ° K ' M ' (29) 

which generalize Eqs. (J23|) - (f21J) above. Note that these results are (independently) perturba- 
tive in each of the coefficients a n and /3 n , because we are assuming that the scalar-field effects 
are negligible with respect to the general relativistic predictions. However, the dominant 
scalar-field corrections may correspond to a large value of index n. 

When solution ()25|) or its approximation (|2Tj) are used to compute the metric coefficients 
a n and j3 n in the above observable predictions, and if we use the huge value of Wq obtained 
from the previous cosmological reconstruction of W(ip), we get the following experimental 
constraint: 

K| < 10- 2xl ° n . (30) 

Now we can speak of fine tuning, and even of hyperfine tuning! This constraint simply means 
that the present value of the scalar field must be exactly at the minimum of the coupling 
function W(ip), otherwise solar-system tests are violated. And since the universe is still 
evolving, the scalar field cannot remain so close to the minimum for more than a fraction 
of a second. Therefore, even if we assumed that Wq = strictly to pass solar-system tests, 
this would not be the case a tiny instant later. In conclusion, we managed to rule out the 
scalar-tensor model A(ip) = 1, V(ip) = and W(ip) ^ 0. It cannot describe an accelerating 
expansion of the universe at present and pass solar-system (and binary-pulsar) tests at the 
same time. 

Of course, this result does not rule out any scalar-Gauss-Bonnet coupling. A model with 
three (or even two) free functions A(<p), V(ip) and W(ip) can obviously pass all present 
tests. [For instance, GR plus a cosmological constant simply corresponds to A(ip) = 1, 
V(<p) = A/2 and W(ip) = 0.] But the presence of a non-constant coupling W(<p) can change 
the physics at small scales, notably in the very early universe (Big-Bang) and for later 
clustering properties. 

The fact that W(<p) induces effects at small scales can be understood by a simple dimen- 
sional argument. Since this function multiplies the square of the curvature in action (J21j) . 
it induces corrections proportional to 1/r 7 (and higher orders) to the Newtonian potential 
in 1/r, and thereby generically dominates at small scales. However, we saw above that 
this quick reasoning can be erroneous in some perturbative but highly nonlinear situations. 
Indeed, if Wq takes very large and negative values, the cosine involved in Eq. (|2Tj) shows 
that ip is always of the order of —Wq/Wq, even for small distances r. One can then prove 
that the (very easily satisfied) condition l^o^ol ^ r2 suffices for all scalar-field effects to be 
negligible in the solar system, even if \Wq \ ~ 10 120 . This remark underlines that nonlinear 
effects can drastically change the intuitive behavior, but let us recall that our cosmologi- 
cal reconstruction above predicted a large and positive value for Wq. In that case, we did 
find that the scalar-Gauss-Bonnet coupling induces large effects at small scales, and even 
exponentially larger than the linear results (|23|) -(|24j ) . 



VI. CONCLUSIONS 



Scalar-tensor theories of gravity are the best motivated alternatives to general relativity. 
Three classes of experimental data give qualitatively different constraints on them. Solar- 
system tests strongly constrain the first derivative of the matter-scalar coupling function 
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A(ip) (i.e., the linear matter-scalar coupling strength ato). Binary-pulsar data forbid large 
and negative values of its second derivative (3q (quadratic matter-scalar-scalar coupling). The 
knowledge of the two cosmological functions Dl(z) and S m (z) suffices to reconstruct the full 
shape of both A(tp) and the potential V(ip) on a finite interval of tp. The knowledge of the 
luminosity distance Dl(z) alone over a wide redshift interval strongly constrains the theories 
if one takes into account solar-system (and binary-pulsar) data, the positivity of the graviton 
and scalar energies, and the stability and naturalness of the models. Future data, provided 
by experiments like the SNAP satellite, will notably allow us to discriminate between GR 
plus a cosmological constant and a potential-free scalar-tensor theory. The possible coupling 
W(ip) of the scalar field to the Gauss-Bonnet topological invariant can be constrained only if 
one takes into account cosmological and solar-system data together. The predictions of the 
model at small distances can depend on highly nonlinear corrections. Of course, a model 
including all three functions A(f), V((p) and W(<p) is experimentally allowed, since GR 
plus a cosmological constant is a particular case. The presence of a scalar-Gauss- Bonnet 
coupling W((p) will generically change the behavior of the theory at small scales (clustering, 
Big Bang). More general tensor-multi-scalar theories, including couplings of the scalars 
to matter and the Gauss-Bonnet invariant, and including non-quadratic "k-essence" kinetic 
terms, can lead to a much richer phenomenology, but at the price of a less natural Lagrangian. 
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